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Conductivity of strongly correlated bosons in optical lattices 
in an Abelian synthetic magnetic field 

A. S. Sajna, T. P. Polak, and R. Micnas 

Solid State Theory Division, Faculty of Physics, A. Mickiewicz University, ul. Umultowska 85, 61-614 Poznan, Poland 

Topological phase engineering of neutral bosons loaded in an optical lattice opens a new window 
for manipulating of transport phenomena in such systems. Exploiting the Bose Hubbard model 
and using the magnetic Kubo formula proposed in this paper we show that the optical conductivity 
abruptly changes for different flux densities in the Mott phase. Especially, when the frequency of the 
applied field corresponds to the on-site boson interaction energy, we observe insulator or metallic 
behavior for a given Hofstadter spectrum. We also prove, that for different synthetic magnetic field 
configurations, the critical conductivity at the tip of the lobe is non-universal and depends on the 
energy minima of the spectrum. In the case of 1/2 and 1/3 flux per plaquette, our results are in 
good agreement with those of the previous Monte Carlo (MC) study. Moreover, we show that for 
half magnetic-flux through the cell the critical conductivity suddenly changes in the presence of a 
superlattice potential with uniaxial periodicity. 

PACS numbers: 03.75.Lm, 05.30.Jp, 03.75.Nt 


I. INTRODUCTION 

The Bose Hubbard model (BHM) is commonly used 
to describe many interesting physical systems e.g. su¬ 
perconductors with short coherence length [Ij, Joseph- 
son junction arrays @1 S or quantum phase transitions 
in cold quantum gases [J, (sj on which the idea of quan¬ 
tum simulations can be realized Recently, attention 
has been paid to the BHM behavior in a strong synthetic 
magnetic field (Ml as well as to its transport prop¬ 
erties [isl - f^ . This research has opened the possibility 
of highly controllable BHM dynamics with explicitly de¬ 
signed kinetics which is a subject of our study. 

Despite of the Josephson junction arrays [H HI up 
to now a strong magnetic field regime has been avail¬ 
able by simulating a vector potential imposed on many 
body wave function of ultra-cold neutral gases. It has 
been realized by engineering the adequate phase of atoms 
when they change their quantum states upon hopping 
through the lattice sites, using e.g. Raman- and photo- 
assisted tunneling (23 - [?7i. sh aking of the lattic e fll . . 
In particular, staggered |25l| and uniform [ 2 ^, l27l| mag¬ 
netic held have been created. Such synthetic magnetic 
held has also been proposed to be generated by combin¬ 
ing quadrupolar potential and modulation of tunneling in 
time [ 2 ^, although it has not been realized yet (see also 
M)- Moreover, the relevant gauge degrees of freedom 
can be precisely experimentally verihable and can cause 
interesting effects like hnite momentum condensate [l3| . 
As follows from the above, the possibility of simulating a 
vector potential provides many opportunities which has 
made it a rapidly expanding area of research. The inter¬ 
est in the simulation has been growing in view of possible 
future application of ultra-cold quantum gases in topo¬ 
logical quantum computation or spintronics (sil. 

If we consider the BHM equipped with orbital effects, 
its complexity considerably increases, as it develops a 
multi-subband structure dependent on the quantity of 


hux per plaquette. To the best of our knowledge the 
conductivity in a strong magnetic held has been consid¬ 
ered mainly numerically in just a few papers. One of the 
main problems in carrying out the related calculations is 
the complex hopping term (Peierls factor) which occurs 
in the BHM Hamiltonian. This is the reason why the 
optical conductivity (OC) has been rarely studied. Y. 
Nishiyama [s^ has analyzed OC in the hard-core limit of 
the BHM with disorder. OC in Josephson junction arrays 
has been also studied within Landau levels framework 
without commensurability effects of magnetic held sm. 
The critical conductivity could be much easier available 
in numerical calculations thanks to the correspondence 
between the BHM at integer hlling and the XY model 
im in which the cases / = 1/2 and / = 1/3 hux per 
plaquette has been studied (^. Using this correspondence 
E. Granato and J. Kosterlitz derived analytically the 
value of critical conductivity for half hux per plaquette. 
Also Y. Nishiyama by applying the exact diagonal- 
ization method has shown that the critical conductivity 
subjected to magnetic held is non-universal. Understand¬ 
ing of such a non-universal behavior gives better insight 
into physics behind the superconductor-insulator phase 
transition phenomena. 

The engineering of the conductivity can be also im¬ 
proved by changing experimentally the value of poten¬ 
tial on adjacent sites. This solution has been recently of 
growing interest because it allows a simple experimen¬ 
tal realizatio n 1^.1^ and generates interesting physical 
phenomena [Hir In particular, tuning the uniaxi- 
ally periodic potential with a magnetic-hux quantum per 
unit cell results in a new possibility of merging cones sim¬ 
ulation in which semi-Dirac point in the Hofstadter spec¬ 
trum emerges [s^ . Recently this possibility has been also 
exploited in the context of time-of-light patterns with and 
without synthetic magnetic field |l3| . 

In this work we present a theory of conductivity valid 
in a strong magnetic held in which only intra-Hofstadter- 
band transitions are taken into account. In the Mott 
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phase for the square lattice, we describe the optical be¬ 
havior using two exemplary values of uniform magnetic 
field: / = 1/2 and / = 1/4 as well as the case with uniax- 
ially staggered potential. We highlight the fact that the 
calculation could be straightforwardly extended to an ar¬ 
bitrary amplitude of the magnetic field and also gauge de¬ 
gree of freedom. For verification of the results obtained, 
a combination of two types of currently available exper¬ 
iments is proposed. Our method is tested for the the 
critical value of conductivity on the Mott insulator - su¬ 
perfluid phase boundary in two dimensions (2D). The 
proposed analytical approach correctly reproduces criti¬ 
cal conductivity in the absence of magnetic field / = 0 
(first calculated in P), / = 1 /2 (MC numerical solution 
in Ref. (2| and analytic one l4(il| but using a model 
correspondin g t o BHM), / = 1/3 (only numerical solu¬ 
tion in Ref. [2| based on MC). We have also extended 
the calculations of critical conductivity over the range of 
arbitrary f = p/q which has not been hitherto reported 
literature and confirmed its non-universal behavior upon 
changes in the amplitude of magnetic field. The influ¬ 
ence of uniaxially staggered potential is also discussed. 
We compare our results with presently available numeri¬ 
cal and experimental data. 

The paper is organized as follows. The BHM in the 
strong uniform magnetic field is reviewed in Sec. IH] In 
Sec. m] we apply our method to study the optical con¬ 
ductivity and its critical value for different strengths of 
synthetic magnetic fields, taking into consideration the 
effects of uniaxially staggered potential. In the last sec¬ 
tion we give a short summary of our results. 


II. MODEL 

The BHM Hamiltonian in standard notation is given 
by 


H = -Y^ Jy (blbj + H.c.) -1- ^ ^ hi (fii - l)-/r ^ fii 

(ij) i i 

where bi (&|) is an operator which annihilate (create) 
boson on site i and hi = b\bi. U and p are on-site boson 
interaction and chemical potential, respectively. Hopping 
integral is denoted by Jij and has a form 

4 = ( 2 ) 

with non-zero isotropic factor Jij = J in respect to ad¬ 
jacent sites. Magnetic field B is introduced by a vec¬ 
tor potential Aq = H(0,x, 0) which is taken in the 
Landau gauge. Further in our calculation we define 
Ba^e*/hc = 27r/ where / = Ba^e*/hc = p/q is a flux 
per plaquette (p and q are coprime integers). / depends 
on the charge of the boson e*, lattice spacing a, Planck 
constant h and speed of light c. It is important to stress 


that in an optical lattice the quantities like B and e* are 
effectively created through tunability of the p/q ratio. 

Using the coherent state path integral for the BHM 
hamiltonian (Eq. (HJ), the partition function can be writ¬ 
ten as follows 

Z[Ao] = y (3) 


5'o 


^ / dr\b*{T)h-^h{T) (4) 

2 0 ^ 

+ ^K{'r)K{'r)h{T)h{T) - ph*{T)b,{T)^ , 


Si [Aq] 


nhP . . 

^ / dr J^jC^ {bUT)bj{T) + c.c.), 

(b> 


( 5 ) 

where the integrals in Eq. (|3]) and Eq. m are taken over 
imaginary time r and /3 = l/ksT. 

To take into account magnetic field in the strong cou¬ 
pling limit of BHM we follow the same procedure as 
in Refs. [H m i.e. we perform the double Hubbard- 
Stratonovich transformation together with cumulant ex¬ 
pansion. In the following, we focus on the Mott phase 
and we approximate effective action to second order in 
the fields 


geff ^ [BlJ [-hGo\^uJn)l + J^(k)] HL, (6) 

kn 


where I is an identity matrix and 

^kn (^k,n5 ^kH-p,n; ; ^k-t-(<?—l)p,n) j ("f) 

with p = (27r/, 0) and uJn = 27rn/tiP is the Matsub- 
ara frequency (n is an integer number). In Eq. ([6]) 
the summation is performed over wave vectors k within 
the first magnetic Brillouin zone where |fcx| ^ n/qa and 
|fcy| ^ 7r/a. GijiiuJn) is the on-site Green function (i.e. is 
local where J = 0) and at the zero temperature limit it 
has the form 


-GniiuJn) = TT — 

a laujri 


no + l 


no 


■ Eno - E, 


rao-l-l 


itUjJr, 


Euq — 1 Ejiq 


( 8 ) 


with on-site energy En^ = —pno + Unoijio — l)/2. no 
is the integer number obtained from the on-site energy 
minimization for a given chemical potential p. The last 
object in ([6|) that requires explanation is J'^(k). It repre¬ 
sents a non-diagonal qxq matrix similar to that obtained 
earlier (4^ 


J^(k)/J = 


Mo 

y^ — ik^U 


_ 

Ml 


^—iknjCL~ 


r^ik-u 


Mg-2 — 

_^-tkya Mq-i 

( 9 ) 
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where = —2cos{kxa + 27raf). The set of eigenvalues 
of the matrix ([5]) leads to the Hofstadter spectrum (d^ . 

It is important to comment here on the notation used 
above. Although b, b* fields introduced in Eq. ([51) have 
identical notation as in Eq. ([3|) they are not the same. 
Namely, those fields in which the quadratic effective ac¬ 
tion Eq. m is evaluated, were introduced during the 
second Hubbard-Stratonovich transformation mentioned 
above. However, both kind of fields have the same corre¬ 
lation function |4l| and for this reason we treat them on 
an equal footing. 

For the effective action in Eq. (I5|) we can find a unitary 
matrix ^^(k) Q that diagonalizes it, i.e. 

geff = -J2{BLy , (10) 

kn 

where = t/,(k) [hGo\iujn)l - J^(k)] t/t(k) 

is an inverse of diagonal Green function and 
= C/q(k)H^^. As shown in Sec. |HT] this diago¬ 
nal form of t/'^(k, has q-bands which are labeled by 
a number. 

The above denotations will be useful for further study. 
In order to simplify the calculations we also set the lattice 
constant a and the reduced Planck constant H to unity. 


Refs. giSI) 

0-4“ (*w) = - 




dre^' 


(52 In Z [A'] 


(5Af (t)(5A“(0) 


A=0 


( 12 ) 

but in our case the conductivity depends upon vector po¬ 
tential Aq and consequently on magnetic field and in this 
dependence the commensurability effects are included. 

Performing the calculations in Eq. (HID we first go to 
the wave vector representation which allows proper incor¬ 
poration of the perturbing vector potential A, through 
the substitution kx ^ kx + j-A. (we assume that A does 
not depend on position) and then we evaluate the deriva¬ 
tives in Eq. HID getting 


0-4“(*w) = {eliJ + -Bxx{i<^) , (13) 

OJ OJ 

where functional averages are taken with the partition 
function from Eq. ([31. ^^nd T\xx{^) are a kinetic and 
current-current correlation function, respectively (we set 
c = 1). With the effective action from Eq. (fTUl) we can 
evaluate Eq. (fT3l) to the form 

Bxx{^^) = (0)), (14) 

^ kk' -lo 


III. CONDUCTIVITY IN MAGNETIC FIELD 


To construct a correct theory of optical conductivity 
in the strong magnetic field described by the BHM, we 
should go beyond the linear response regime. This situ¬ 
ation is more c omp licated than that previously studied 
(e.g., see Ref. [^) due to the presence of a synthetic 
magnetic field Aq which significantly modifies the single¬ 
particle spectrum. This modification implies that the 
commensurability field effects become important. As a 
consequence, the amplitude of hopping term is a complex 
number (Eq. ([SD) and the boson field is described by 
more than one component (Eq. (0). To overcome these 
difficulties in calculation of the transport properties we 
assume that such an initial state will be subtly affected 
by an additional field | A| ^ l-A-ol, which is responsible for 
generation of optical conductivity in the linear response 
regime. Therefore, if we want to study the response of 
the system to a strong magnetic field, we simply add to 
the Aq some small perturbation of this field in the form 
of a vector potential A i.e. 



Aq • d\ —^ 



Aq ■ dl -f 



A-dl 


( 11 ) 


and use the linear response theory as a starting point 
with respect to the quantity A. This leads to the well 
known expression for the optical conductivity (e.g., see 


^kin 


1 

N 


9-1 

K+ap{0)bk+cp{0), 

k a=0 


(15) 


9-1 


jg(r) = e* ^ [9fc^e“(k;p)] &k+cp(i‘)^k+ap(T), (16) 


q;=0 


where e“(k;p) is an eigenvalue of Jg(k) matrix defined 
in Eq. (0 and the current is given by Eq. (ITSl) . Eq. 
(fT31) is the Kubo formula valid in strong magnetic field 
(MKF) in which only intra-Hofstadter-band transitions 
are taken into account. This method, in contrast to that 
presented in Ref. @ (see also Refs. [13, EE E3“E3)> 
takes into account the commensurability effects of the 
magnetic field and covers the whole range of k and uj 
dependence. 

For the simplest case, if we take the usual square lat¬ 
tice dispersion relation in the absence of magnetic field 
e“(k;0) = — 2J {coskx -f cosky), one gets from Eq. (fTSl) 


= ^^-^^^‘^coskx{bl{0)bk{0)} (17) 


E / 


4(e*)^j2 1 

UJ N - Jr) 

k.k' ^ 

sin(fcj,)sin(fc^)(5k(T)6k(r)6k,(0)6k'(0)) 

which recovers a well known result Ejl. 
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w/U 



rewrite Eq p3ll to the form 




9-1 


= - 


w PN 
{e*f 1 




kn q;=0 
g-1 

X 

kn CK^O 
2 


EE 


( 18 ) 


w /3A1 

C*-\J 

[5fcx e, (k;p)] ^ Gta (k, iuJn)Gac (k, *Wn + iuj). 
Mott insulator Green’s function is 
Go(*Wri) 




1 - e“(k;p)Go(ia;„) 


(19) 


<(k;p) 


1 - z^{k;p) 


iWn - E,+(k; p) iuJn - Eq (k; p) ’ 

where the weight and dispersion of quasiparticles have 
the form 


<(k;p) = 


E“+(k;p) + p + t/ 
E^+{k-p)-E^-{k-p) 


( 20 ) 


Eq^[k-p) = -p + G ± ^^9(k;p) 

( 21 ) 

A“(k;p) = ^[e“(k;p)] V4e“(k;p)[/ +1^2 

( 22 ) 

In the following we are interested in the real and reg¬ 
ular part of optical conductivity given by 

27r 

Recr^x (w) = {e*f ^ E E (k;p))^ (23) 

Q—0 k 

X {nB [£’“"(k;p)] - ns [£’“+(k;p)] } 

X [1 - 2;“(k;p)] z’^{k]p)5 (uj"^ - (A“(k;p))^) , 


Figure 1: Real part of optical conductivity in the two- 
dimensional square lattice in Mott phase (first lobe) is 
sketched in ag units (ag = (e*)^ /h is a, quantum conduc¬ 
tance). From the top - / = 0 (a), / = 1/2 (b), / = 1/4 
(c). Conductivity is plotted at zero temperature for different 
value of J/U within the first lobe. 


which has been obtained by evaluating Matsubara sum 
in Eq. (I18II and performing standard analytical continu¬ 
ation iuj ^ uj + iS d^. The symbol 6{x) is a Dirac delta 
function and hb^x) is Bose-Einstein distribution func¬ 
tion. Using density of states and taking zero-temperature 
limit, equation (1231) can be rewritten as follows 


A. Optical conductivity in the BHM 


= 27rVQ^ 5 “[m'*(w);p] ,(24) 

“ s={+.-} 


1. The uniform field 


'^a 


b;p] = Pnivip)- 


J [zq{v]p) - 1] 2:“(u;p) 


,(25) 


U J 4no (no + 1) + (w/U) 


Now we are interested in the OC in the Mott insulator 
phase. Using the fact that within the action m the four 
point correlation function in Eq (1141) is factorized, we can 


“*(“) = 7<2"o + 1) p/'/O',)) <■*) 
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Figure 2: The wave vector k dependence of the quasi-particle 
energy dispersions for different amplitudes of uni¬ 

form magnetic field / and staggered potential A in the first 
magnetic Brillouin zone. Plots a, b, c correspond to the Figure 
[T]a, b, c, respectively. Additionally, plots d and e correspond 
to Fig. [4] and Fig. [5] respectively. This correspondence is 
revealed when choosing the same physical parameters for the 
same plotting style in the relevant figures. 

where the density of states is given by 

= . (27) 

k 

For magnetic fields considered here (i.e. / = 0, 1/2, 1/4) 
the exact calculations of density of states in 2D was pos¬ 
sible in terms of complete elliptic integrals (see Appendix 

ESI- 

In Fig. [T] we show the real part of OC at zero tem¬ 
perature in the Mott phase for different values of syn¬ 
thetic magnetic field. It is worth noting that its be¬ 
havior reflects the tight binding dispersion of the lattice. 
As expected the OC is gradually broadened at the cost 
of vanishing gap when J/U increases. Interestingly, for 
/ = 1/4, the contribution of the lowest frequency peak 
becomes much more significant when the J/U ratio is 
tuned up. 



Figure 3: Optical conductivity for different values of magnetic 
field. The parameters fi/U and J/U are chosen on the phase 
boundary where Mott insulator - superfluid phase transition 
takes a place. The plot shows that the critical value of con¬ 
ductivity for the / = 1/2 (1/4) when a; —>■ 0 is two (four) 
times higher than in the case when magnetic field is absent. 
Rea^x was plotted in ag units. 


For Lo = U,we observe that the existence of conductiv¬ 
ity at this point directly depends on the spectrum weight 
value of the tight binding dispersion e“(k;p). In particu¬ 
lar, if the spectrum weight of e“ (k; p) at the center of the 
band is zero the conductivity also disappears. Such a sit¬ 
uation for w = [/ is satisfied when f = 1/2, 1/4 (even q) 
but for odd q we should observe metallic behavior. This 
special behavior of OC at this point is directly related to 
the Dirac cones appearing in quasi-particle energy dis¬ 
persion A“^(k;p). The corresponding dispersions are 
plotted in Fig. 2 a — c for the relevant set of parame¬ 
ters. In the Sec. IIII BI we suggest an experiment to check 
this conjecture explicitly, because the quantities like the 
on-site interaction strength U and boson hopping am¬ 
plitude J are fully controllable parameters in ultra-cold 
gases loaded on optical lattice. 

Within the above framework we can also probe the 
critical value of conductivity at the tip of the lobe (for 
determination of the critical value oi p/U and J/U see 
Ref. [HI). For this range of parameters the OC is as 
shown in Fig. [31 The critical value of conductivity for 
/ = 1/2 (1/4) when a; —>■ 0 is two (four) times higher 
than the value with magnetic field absent. These results 
are in agreement with those derived in Sec. IIII Cl 


2. The uniform field with staggered potential 

In the following we use the proposed framework of OC 
(see Eq. (fT31) ') to study the transport phenomena under 
the influence of uniaxially staggered potential. The sit¬ 
uations with the synthetic magnetic field absent (/ = 0) 
and with its value described by half flux per plaquette 
(/ = 1/2) are considered. The latter special case is a 
subject of current interest (s^. In both cases of OC i.e. 
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Figure 4: First (second) panel shows yy (xx) component of the real part of the optical conductivity for different values of A 
in the absence of magnetic field. On the density sketch the lighter color indicates higher amplitude of the conductivity (black 
area correspond to the insulator behavior). All zero-temperature plots show the situation where the uniaxially staggered field 
is subjected to the x axis. The ratio of hopping amplitude and on-site interaction energy is J/U = 0.03 and the conductivity 
was plotted in ctq units. 


/ = 0 and f = 1I2 the staggered potential is controlled 
by the parameter A , which is included in the Hamil¬ 
tonian ([T]) by an additional term in the form 

, (28) 

ixiy 

where A = AI2J and ix, iy enumerate the positions of 
lattice sites along the x and y axis, respectively. 

Fig. 0] and show the frequency dependence of OC 
for / = 0 and / = 1/2, respectively. We are interested 
in and component in which the potential 

from site to site is varied in the x direction (to see the 
expressions used in the calculations of OC see Appendix 

lYH). 

The data presented in Fig. 0] and 0] imply a similar 
behavior of OC when A parameter is alternated. For 
example, on the basis of the staggered potential values 


with respect to its value for A = 0, we conclude that 
it has greater impact on the xx component of OC than 
on yy one. Besides the qualitative difference, we also 
observe a smaller amplitude of the OC in direction x 
than in that perpendicular to x within xy plane. This 
behavior could be simply attributed to the variation in 
the potential in this particular direction (i.e. cc). While 
the yy component of the OC does not exhibit any special 
difference along y axis (for chosen strip of sites the A 
is constant). Moreover, we show that the increase in A 
causes broadening a frequency dependence of OC. 

In agreement with the conclusion drawn in Sec. IIII A II 
the insulator behavior of the OC for w = t/ (Fig. 0]) is 
still maintained for / = 1/2. In contrast to the situa¬ 
tion with no magnetic field (see Fig0|) the gap naturally 
arises when the A parameter exceeds one (the single¬ 
particle spectrum also exhibits a similar behavior). This 
gap-like behavior is indeed observed in the quasi-particle 
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energy dispersion presented in Fig. [5]d and e. 

Interestingly, the weight of the OC close to uj = U for 
/ = 1/2 is greater than for / = 0. 

It is worth noting that the complex behavior of the OC 
in the direction of applied uniaxially staggered potential 
gains pronounced peaks. This should be easily observed 
in the EAR experiment (see Sec. IIII Bl) . 


B. Connection to experiment 


Recently, A. Tokuno and T. Giamarchi [l^ have pro¬ 
posed a spectroscopic technique for cold atoms which 
is able to extract current-current correlation function 
(Eq- (HHl)- This function is proportional to OC 
(see Eq. (fT3l) ), which offers a possibility to probe 
the transport phenomena in the thermodynamic limit. 
Using the energy absorption rate techniques (EAR) such 
a goal could be achieved by phase modulation of optical 
lattice. 

Namely, a vector potential Aq could be created in dif¬ 
ferent experimental configurations (2^, [2^, [2^. In our 
work we investigate the Landau gauge Aq = B{0,x,0) 
which generates a uniform magnetic field. Such a uni¬ 
form field has been generated recently in but using 
another type of effective vector potential. On the other 
hand, a small perturbing vector potential A' could be 
generated by phase modulation of optical lattice M- If 
we consider a 2D system with x and y axes, we can gener¬ 
ate a synthetic electric field Ex by modulating the phase 
in the x direction. This situation could be mathemati¬ 
cally inferred from the exchange of a stationary optical 
lattice potential I4p(r) = cos^ (fea;a;)-|-cos^ (^yl/) to a time 
dependent one 14p (r, t) where x ^ x — fx cos iujt) (fx is 
the strength of modulation which should be much smaller 
than a lattice constant). Such a phase modulation (PM) 
could be realized by e.g. a recently proposed phase con¬ 
troller [s^l- This leads us to the expression where EAR 
is given by [l^ 


Rpuito) = -iwVJSn,,(w) , (29) 

in which Ilxx{uj) corresponds to Itxxi^x)) from Eq. o 
but with exchange e* —>■ M (M is effective mass of an 
atom) [l^ . To ensure a linear response regime and no 
dynamical phase transition the condition w/a, -C 1 should 
be satisfied (sil. Is^. The calculation of i?pM(w) is made 
following a procedure similar to that used in the OC case 
where the current-current correlation function Ytxx{i<jj) 
was also considered (see Eq. 0). Fig. ® presents 
a plot of Rpm/^ qfx for tho uniform magnetic field of a 
strength / = 0, 1/2, 1/4 on the square lattice in the 
zero temperature limit (here dq = M'^/h is an effective 
quantum resistance). We see that the factor changes 
significantly the weight of the response in comparison to 
the OC given in Fig. [l]and therefore the higher frequency 
peaks give a grater contribution to the absorbed energy 


rate. The similarity of the shape of current-current cor¬ 
relation function and tight-binding dispersion in a strong 
magnetic field could be an indirect method of checking 
the Hofstadter spectrum in the BHM system [1,|45 where 
the center of the band is located around lu = U. 

The plots of EAR for the OC with staggered potential 
will be analogous to the case of a uniform field (discussed 
above) so we omit its graphical representation. 

Summarizing, the EAR technique is directly related 
to OC and may act as a probe of transport phenomena 
using ultra-cold quantum gases. It is worth pointing out 
that the phase modulation is independent of the strength 
of the lattice potential in contrast to the amplitude mod¬ 
ulation method [T^ . 


C. Critical conductivity at the MI-SF phase 
boundary 


1. The uniform field 


Up to now the problem of critical conductivity in 2D 
on the Mott insulator - superfluid phase boundary in 
strong magnetic field has been rarely studied because 
of its complexity. The amplitude of the hopping term 
with a complex factor (see Eq. ([S])) is the reason why 
in Ref. instead of BHM at integer filling, the frus¬ 
trated XY model was investigated using the MC numer¬ 
ical method. Another approach to the problem of criti¬ 
cal conductivity in magnetic field has been proposed in 
Ref. [U, but the magnetic field considered there was 
weaker (i.e. / = 1/20) than that the field we studied 
and their calculations were performed in the hard-core 
limit (the author used exact diagonalization method to 
overcome the difficulties related to the complex hopping 
term). Within the approach presented in this paper we 
can perform analytical analysis and study critical behav¬ 
ior of conductivity in a much wider range of magnetic 
fields. Namely, we show that for commensurate value of 
/, critical conductivity depends only on the number of 
minima located in the first reduced magnetic Brillouin 
zone. 

To describe the critical conductivity at the tip of the 
Mott lobe El we compute the OC (Eq. (IT^ l close to 
the phase boundary. To do that, in the following calcula¬ 
tions, we consider only the real part of (w) that gives 
a finite frequency contribution, namely the part of OC 
which consists of the current-current correlation function 


ReCT:^“ {ioj) = —Re 


_Li V V X 

Rq h/3N huj h 


(30) 


[9fcxeg(k;p)]^eL(kuWn)C?L(kuw„ + iw)| , 


where Rq = hj = Oq is quantum resistance (for 
Cooper pair e* = 2e and Rq ~ 6,45 /cD) and we restore 
the constant h to introduce it in quantum resistance Rq. 
Eq. (1301) also contains the singular part of OC, but since 
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Figure 5: Analogous physical situation as in Fig. Q but with non-zero amplitude of uniform magnetic field (half magnetic-flux 
per unit cell). 
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we are interested in the Mott phase, we neglect this con¬ 
tribution further on. 

Now, using the effective action from Eq. (HOD, and ap¬ 
plying the Ginzburg-Landau ( GL) like method for calcu¬ 
lation of critical conductivity [44| . we evaluate Eq. 
in order to obtain its dependence on / = p/q parameter. 


Following this procedure, we expand the action Eq. m 
to the second order in frequency using the expression 

1 - e“(k; p)Go{iuin) = flk - buifu^n - Ck (fSw„)^ , (31) 
with 


Ok = 1 - eg(k;p)Go{iuJn = 0) 


(32) 
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6k = e^{k;p)drGQ{r)\r=o , Ck = ie“(k;p)a^Go(r)|r=o • 

Next, in calculating the critical conductivity within the 
GL action, we should assume the proper ground state 
behavior. From all set of g—band energy dispersion, 
we choose the lowest one which correctly reproduces the 
phase transition. This band contains q GL modes in the 
first magnetic Brillouin zone (MBZ) [ij, which allows 
description of the critical behavior of the BHM close to 
the phase boundary. Going further we perform the sum¬ 
mation over Matsubara frequencies and take the limit 
T —>■ 0, which reduces Eq. (1^ to 




1 \ ^ (k — Q) 2 / 2 ^k 


(33) 

where Ak = y/bl + dokCk, We// = 1/ e“ (k; p) | k=Q 
and Q are locations of the minima in MBZ. Glose to 
the phase boundary only the momenta around Q bring 
a contribution to conductivity, therefore if the minimum 
of e“(k;p) is located at k = Q we can simply expand Uk 
from Eq. (15^ to the second order 


O-k 


(k-Q)^ 

2zJmeff 


(34) 


where z = e“(Q;p) for chosen p, q. In further calcula¬ 
tions we assume that i9^^e“(k;p) = 9^^e“(k;p). Einally, 
for the point in the phase diagram which is close to the 
tip of the lobe, we get 


8Rq ^ U}^ 

(35) 

where 0(a;) is a step function being non-zero for x ^ 0. 
The above expression describes the behavior of opti¬ 
cal conductivity close to the phase transition. It has 
non-vanishing amplitude when the applied frequency is 
equal to 6q/cq or is higher than this value. Denot¬ 
ing a* = naqlS = ir/SRq and considering the tip of 
the lobe, where 6q = 0 is the critical conductivity, 
Re(f;^“(a;) = CTc,/ takes the simple form 

CTc./ = qcr* . (36) 

It is important to notice that theory presented here is 
valid for the second order phase transition, e.g. this con¬ 
dition is satisfied for the cases / = I/2, 1/3 [3- 

To discuss the result from Eq. (l36l) we firstly recall that 
for the simplest case g = 1 where there is no magnetic 
field, the result for g = 1 confirms the results presented in 
Refs. ( 1 , 111 , 13 . With p/q = 1/2 we have acp /2 = 2cr*. 
This result agrees with the analytical solution given in 
Refs. [ 11 , 1411 , where the XY model was used. Also in the 
MG study for / = 1/2 dl a value of critical conductivity 
is 1,82 times higher than for / = 0 [111. Results of an 



f 

normalized critical conductivity ac,f/o'c,o 


MKF (here) XY model 

MC JJA 

experiment 

1/2 

2 2 

1,82 « 2 

1/3 

3 

2,91 « 3 

1/q 

q 

« q 


Table I: Comparison of critical conductivity for different mag¬ 
netic field configuration calculations 0, IH, IH, Ull and exper¬ 
imental measurements [H, (Hi . 


exp eriment conducted in Josephson junction arrays [H . 
[Hi also show a similar behavior. If we consider / = 1/3, 
it qualitatively agrees with the Monte Garlo result in Ref. 
d, where Gha and Grivin obtained 2, 91 higher value of 
critical conductivity than in the absence of a magnetic 
field. The authors of the experimental work, in Refs. 
[H, nil also discuss such a scenario. It is worth adding 
that the authors of Ref. [Hi have speculated about a 
similar result for cFcp/s i-e. <^ 0 , 1/3 = 3(Tc,o- Namely, they 
have suggested that at least for low order rationals of 
p/q critical conductivity could satisfy Eq. (1551) but they 
have carried out explicit calculation only for f = 1/2. In 
contrast, we showed this behavior by analytical methods 
for arbitrary / within a clear mathematical framework. 
The above considerations are summarized in Table ID It 
is worth mentioning that in three dimensions a trivial 
solution CTc,/ = 0 is obtained, known before only for the 
case of / = 0 (e.g, see Hi)- 

Eor Josephson junction arrays (JJA) it seems that the 
critical conductivity is proportional to q/p (H but this 
linear behavior is inferred from a small number of ex¬ 
perimental points with a large error margin. Therefore, 
such a dependence is still an open question. Moreover, 
in (JJA) we should take into account that arrays are not 
perfect and their parameters differ through the network. 
Also, the measurements are performed at finite temper¬ 
atures. For example, if we consider a disorder we should 
expect that this effect suppresses the value of critical con¬ 
ductivity BlH- It is worth mentioning here that also 
long range interactions which we neglected could have a 
significant impact [H- Besides, it seems that CTcJ should 
depend on p. Hence our results (1361) within the approxi¬ 
mations used in this paper should be at least appropriate 
for p Ri I. 


2. The uniform field (f=l/2) with uniaxially staggered 
potential 

To show the importance of translation symmetry 
breaking by the uniaxially staggered potential in a spe¬ 
cial case / = 1/2 [H j we analyze the critical behavior of 
conductivity at the phase boundary. 

Following the same procedure for the critical conduc¬ 
tivity which led us to Eq. (1361) we simply observe that for 
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I I I I I I 

"D" 

I I I I I I 


Figure 7: Above we sketch density plot of the lowest band 
energy in the Mott phase (white rectangle describe first 
magnetic Brillouin zone (MBZ)). The range from darker to 
brighter color is assigned to the lowest and highest value of en¬ 
ergy spectrum respectively. Pictures have A = 0 and A = 0.7 
(from left). We see explicitly disappearance one of the two 
minima (black color) on the center of the MBZ. 

A = 0 (for definition of parameter A see Sec. IIIIA2I) we 
get acj = 2a*. But for A > 0, the situation is changed 
significantly. Analysis of the spectrum of quasi-particles 
in the Mott phase reveals that one of the two minima 
disappears in the first magnetic Brillouin zone for the 
non-zero value of A (see Fig. [21). Therefore, there exists 
only the one lowest energy Ginzburg-Landau mode which 
effectively recovers the critical conductivity as when there 
is no magnetic held. Hence, such an abrupt change in the 
critical value of CTc,/ from 2tT* to cr*when the staggered 
potential is turned on, could be an interesting effects on 
its own right. In addition, no variation in the critical 
conductivity for / = 0 was observed. 

The best options to study the presented effects in ex¬ 
periments that could be made nowadays is the use of 
optical lattices with ultra-cold atoms whose high con¬ 
trollability provides a better road to connect experiment 
and theory. The continuous progress in ultra-cold quan¬ 
tum gases in the near future also will allow getting into 
the temperature regime where high precision measure¬ 
ment of critical conductivity will be possible. This can 
verify the above results and allow omission of additional 
effects which occur in standard solid state devices like 
Josephson junction arrays. The possibility of such mea¬ 
surements has been very recently discussed in Ref. (^ . 

The analysis made in Sec. IIII Cl permits a better un¬ 
derstanding of the superconductor-insulator phase tran¬ 
sition mechanism. To illustrate this, we shortly explain 
below the fact that the critical conductivity should be 
affected by the applied magnetic held. Namely, we know 
that the scenario of critical resistance (i.e. I/ctc,/) is as¬ 
signed to the vortex and boson flow ing through the sys¬ 
tem at the critical transition point [46|. Such a descrip¬ 
tion is possible from the duality transformation 
(e.g. for superconductor-insulator transition induced in 
a fermionic system these bosons are Cooper pairs with 
short coherence length Hi). Therefore, application of a 
magnetic held at least affects the behavior of vortices. 


which changes the critical resistance. In our calculations 
the magnetic held is effectively incorporated into the the¬ 
ory through the tight binding dispersion relation e“ (k; p) 
which gives g—bands spectrum with q minima in the low¬ 
est energy level, what hnally changes CTc,/- If we addi¬ 
tionally consider the staggered potential, an analogous 
prediction could be made. 

IV. SUMMARY 

The analysis of conductivity in BHM in a strong mag¬ 
netic held is challenging problem due to complex hopping 
term. In particular, up to now its optical dependence was 
out of reach in Monte Carlo study. Therefore our theory 
expands the area in which numerical methods are used. 

Namely, we have proposed the magnetic Kubo formula 
which is valid for an arbitrary hux pattern where com- 
mensurability effects of a magnetic held are included. 
Within this framework, we have calculated the opti¬ 
cal conductivity in the Mott phase of the Bose Hub¬ 
bard model and considered its critical value. To check 
our results we have proposed to compare them with 
presently available experiment in ultra-cold quantum 
gases in which current-current correlation function in a 
uniform magnetic held could be probed. Such a con¬ 
nection of experiments and theory could open a new av¬ 
enue to study transport phenomena in highly controllable 
magnetic held where geometry of the lattice can be easily 
manipulated. Moreover, for the case of critical conduc¬ 
tivity we have shown its dependence on the topology of 
the single-particle spectrum and obtained solution which 
are in good agreement with presently available numerical 
and experimental data. 

The method presented here can be extended over 
many-body systems in which the strong magnetic held 
plays a signihcant role. 
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V. APPENDIX 

Here we present the quasi-particle energy spectrum 
and densities of states for conductivity (DOSc) in 2D 
described in Sec. MM We use K. and £ to denote the 
hrst and second complete elliptic integral, respectively. 

A. Uniform magnetic field for a two-dimensional 
square lattice 
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1. f=0 


The dispersion relation for the square lattice is 

e“(k; 0) = —2J (cos+ cos ky ), 

where the DOSc (Eq. (l27ll is given by (55l | 


Pi{E;0) 




(37) 


(38) 


a. f=l/2 


The dispersion relation for / = 1/2 is built from two sub-bands ±2 y^cos^ kx + cos ky and consequently the 

DOSc is 


P2{E-,1) 




where 0^ is non-zero step function within each q-band. 


(39) 


3. f=lU 


The form of dispersion relation / = 1/4 is expressed by four sub-bands ±y^d — -^12-1-2 cos {4:kx) + 2 cos (4^^ 
and ±y^4-|- -^12-1-2 cos (dfe^^) -|- 2 cos (4fcy), then we get 


p2iE;i) = 


40c 


7t‘^\E^-4\\E\ 


' 1 - 


A-8E^ + E^ 


4 - 8^2 -h E‘^ 


K. 


' 1 - 


A-8E^ + E^ 


iM 


B. Uniform magnetic field for a two-dimensional 
square lattice with uniaxially staggered potential 


1. f=0 


a/rp 1'V 

P2iE]l} = — 

—r^ 


20a 

Jo 


dx A -- A . (42) 


\/l — x'^ 


The form of tight binding dispersion has two sub-bands 
( 3 ^ —2Jcosky ± 2JVcos^ kx + and the appropriate 
DOSc is given by 

- for XX component of optical conductivity a^°{uj): 


20 

p^iE;l) = ^ 


20a 

Jo 


dx ■ 




A2 


^l-d 


(41) 


for yy component of optical conductivity cr^°(a;): 


a. f=l/2 

The form of tight binding dispersion has two sub¬ 
bands [s^ ±2'y/cos2 kx + (cos ky — A)^ and the appro¬ 
priate DOSc is given by 

- for XX component of optical conductivity 


P^{E;1) = 


40„ 


>^\E\ 


-1-A 


1 - (x + A)^ 


(43) 
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- for yy component of optical conductivity 


P2iE;i) 
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